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Abstract 

Many nonparametric regressors were recently shown to converge at rates that de- 
pend only on the intrinsic dimension of data. These regressors thus escape the 
curse of dimension when high-dimensional data has low intrinsic dimension (e.g. 
a manifold). We show that fc-NN regression is also adaptive to intrinsic dimen- 
sion. In particular our rates are local to a query x and depend only on the way 
masses of balls centered at x vary with radius. 

Furthermore, we show a simple way to choose k = k(x) locally at any x so as to 
nearly achieve the minimax rate at x in terms of the unknown intrinsic dimension 
in the vicinity of x. We also establish that the minimax rate does not depend on a 
particular choice of metric space or distribution, but rather that this minimax rate 
holds for any metric space and doubling measure. 



1 Introduction 

We derive new rates of convergence in terms of dimension for the popular approach of Nearest 
Neighbor (fc-NN) regression. Our motivation is that, for good performance, fc-NN regression can 
require a number of samples exponential in the dimension of the input space X. This is the so-called 
"curse of dimension". Formally stated, the curse of dimension is the fact that, for any nonparametric 
regressor there exists a distribution in M. D such that, given a training size n, the regressor converges 
at a rate no better than n,~ 1 /°(- D ) (see e.g. fl~l|Z)). 

Fortunately it often occurs that high-dimensional data has low intrinsic dimension: typical examples 
are data lying near low-dimensional manifolds (3] |4] [5] . We would hope that in these cases non- 
parametric regressors can escape the curse of dimension, i.e. their performance should only depend 
on the intrinsic dimension of the data, appropriately formalized. In other words, if the data in R D 
has intrinsic dimension d « D, we would hope for a better convergence rate of the form nT 1 / ^ 
instead of n -1 /°( D ) . This has recently been shown to indeed be the case for methods such as kernel 
regression [6], tree-based regression Q and variants of these methods [8|. In the case of fc-NN 
regression however, it is only known that 1-NN regression (where fc = 1) converges at a rate that de- 
pends on intrinsic dimension |9). Unfortunately 1-NN regression is not consistent. For consistency, 
it is well known that we need fc to grow as a function of the sample size n iflOl . 

Our contributions are the following. We assume throughout that the target function / is Lipschitz. 
First we show that, for a wide range of values of fc ensuring consistency, fc-NN regression converges 
at a rate that only depends on the intrinsic dimension in a neighborhood of a query x. Our local 
notion of dimension in a neighborhood of a point x relies on the well-studied notion of doubling 
measure (see Section |231 ). In particular our dimension quantifies how the mass of balls vary with 
radius, and this captures standard examples of data with low intrinsic dimension. Our second, and 
perhaps most important contribution, is a simple procedure for choosing fc = k(x) so as to nearly 
achieve the minimax rate of O (n~ 2 / ( 2+d )) in terms of the unknown dimension d in a neighborhood 
of x. Our final contribution is in showing that this minimax rate holds for any metric space and 
doubling measure. In other words the hardness of the regression problem is not tied to a particular 
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choice of metric space X or doubling measure p, but depends only on how the doubling measure p 
expands on a metric space X. Thus, for any marginal p on X with expansion constant 9 (2 rf ), the 
minimax rate for the measure space (X, p) is D, {n~ 2 ' ^ 2+d '\ 

1.1 Discussion 

It is desirable to express regression rates in terms of a local notion of dimension rather than a global 
one because the complexity of data can vary considerably over regions of space. Consider for ex- 
ample a dataset made up of a collection of manifolds of various dimensions. The global complexity 
is necessarily of a worst case nature, i.e. is affected by the most complex regions of the space while 
we might happen to query x from a less complex region. Worse, it can be the case that the data 
is not complex locally anywhere, but globally the data is more complex. A simple example of this 
is a so-called space filling curve where a low-dimensional manifold curves enough that globally it 
seems to fill up space. We will see that the global complexity does not affect the behavior of fc-NN 
regression, provided fc/n is sufficiently small. The behavior of fc-NN regression is rather controlled 
by the often smaller local dimension in a neighborhood B(x, r) of x, where the neighborhood size 
r shrinks with fc /n. 

Given such a neighborhood B(x, r) of x, how does one choose fc = k(x) optimally relative to the 
unknown local dimension in B(x, r)l This is nontrivial as standard methods of (global) parameter 
selection do not easily apply. For instance, it is unclear how to choose fc by cross-validation over 
possible settings: we do not know reliable surrogates for the true errors at x of the various estimators 
{fn,k{x)} k ^ n y Another possibility is to estimate the dimension of the data in the vicinity of x, and 
use this estimate to set fc. However, for optimal rates, we have to estimate the dimension exactly and 
we know of no finite sample result that guarantees the exact estimate of intrinsic dimension. Our 
method consists of finding a value of fc that balances quantities which control estimator variance and 
bias at x, namely 1/fc and distances to x's fc nearest neighbors. The method guarantees, uniformly 
over all x € X, a near optimal rate of O (n~ 2 /( 2+d ') where d = d(x) is exactly the unknown local 
dimension on a neighborhood B(x,r) of x, where r —> as n —> oo. 

2 Setup 

We are given n i.i.d samples (X, Y) = {(Xi, from some unknown distribution where the 

input variable X belongs to a metric space (X, p), and the output Y is a real number. We assume 
that the class B of balls on (X, p) has finite VC dimension Vg. This is true for instance for any 
subset A" of a Euclidean space, e.g. the low-dimensional spaces discussed in Section [231 The VC 
assumption is however irrelevant to the minimax result of Theorem[3] 

We denote the marginal distribution on X by p and the empirical distribution on X by p n . 
2.1 Regression function and noise 

The regression function f(x) = E [Y\X = x] is assumed to be A-Lipschitz, i.e. there exists A > 
such that Vx, x' e X, \ f (x) - f(x')\ < \p(x,x'). 

We assume a simple but general noise model: the distributions of the noise at points x € X have 
uniformly bounded tails and variance. In particular, Y is allowed to be unbounded. Formally: 

V«5 > there exists t > such that sup P Y \x=x (\Y - f(x)\ > t) < S. 

We denote by ty(S) the infimum over all such t. Also, we assume that the variance of (Y\X = x) 
is upper-bounded by a constant a Y uniformly over all x € X. 

To illustrate our noise assumptions, consider for instance the standard assumption of bounded noise, 
i.e. \Y — f(x)\ is uniformly bounded by some M > 0; then V<5 > 0, ty{5) < M, and can thus be 
replaced by M in all our results. Another standard assumption is that where the noise distribution 
has exponentially decreasing tail; in this case V<5 > 0, ty(S) < 0(ln 1/(5). As a last example, in the 
case of Gaussian (or sub-Gaussian) noise, it's not hard to see that V<5 > 0, ty(S) < 0(y/\nl/S). 



2 



2.2 Weighted fc-NN regression estimate 

We assume a kernel function K : R + h-> K + , non- increasing, such that K(l) > 0, and K(p) = 
for p > 1. For x € A", let rfe, n (a;) denote the distance to its fc'th nearest neighbor in the sample X. 
The regression estimate at x given the n-sample (X, Y) is then defined as 



2.3 Notion of dimension 

We start with the following definition of doubling measure which will lead to the notion of local 
dimension used in this work. We stay informal in developing the motivation and refer the reader to 
iTTTl H"2l rm for thorough overviews of the topic of metric space dimension and doubling measures. 

Definition 1. The marginal p is a doubling measure if there exist Cdb > such that for any x £ X 
and r > 0, we have p(B{x,r)) < Cdbp(B(x, r/2)). The quantity Cdb is called an expansion 
constant of p. 

An equivalent definition is that, p is doubling if there exist C and d such that for any x £ X, for any 
r > and any < e < 1, we have p(B(x, r)) < Ce~ d p(B(x, er)). Here d acts as a dimension. It 
is not hard to show that d can be chosen as log 2 Cdb and C as Cdb (see e.g. 03)). 

A simple example of a doubling measure is the Lebesgue volume in the Euclidean space R rf . For any 
x £ M. d andr > 0, vol(B(x,r)) = vol (B(x, 1)) r d . Thus vol (B(x, r)) / vol (B(x, er)) = e~ d 
for any x £ M. d , r > and < e < 1. Building upon the doubling behavior of volumes in R d , 
we can construct various examples of doubling probability measures. The following ingredients are 
sufficient. Let X c M. D be a subset of a d-dimensional hyperplane, and let X satisfy for all balls 
B(x, r) with x £ X, vol (B(x, r) fl X) = Q(r d ), the volume being with respect to the containing 
hyperplane. Now let p be approximately uniform, that is p satisfies for all such balls B(x,r), 
p(B(x, r) n X) = 6(vol (B(x, r) n X)). We then have p(B{x,r))/p(B(x,er)) = 0(e- d ). 

Unfortunately a global notion of dimension such as the above definition of d is rather restrictive as 
it requires the same complexity globally and locally. However a data space can be complex globally 
and have small complexity locally. Consider for instance a d-dimensional submanifold X of M D , 
and let p have an upper and lower bounded density on X. The manifold might be globally complex 
but the restriction of p to a ball B(x, r), x £ X, is doubling with local dimension d, provided r is 
sufficiently small and certain conditions on curvature hold. This is because, under such conditions 
(see e.g. the Bishop-Gromov theorem fl3)), the volume (in X) of B(x, r) n X is <d(r d ). 

The above example motivates the following definition of local dimension d. 

Definition 2. Fix x £ X, and r > 0. Let C > 1 and d> 1. The marginal p is (C, d) -homogeneous 

on B(x, r) if we have p(B(x 7 r')) < Ce~ d p(B(x, er')) for all r' < r and < e < 1. 

The above definition covers cases other than manifolds. In particular, another space with small local 
dimension is a sparse data space X C M. D where each vector x has at most d non-zero coordinates, 
i.e. A" is a collection of finitely many hyperplanes of dimension at most d. More generally suppose 
the data distribution p is a mixture t^iPi of finitely many distributions pi with potentially differ- 
ent low-dimensional supports. Then if all pi supported on a ball B are (C, d) -homogeneous on B, 
i.e. all have local dimension d on B, then p is also (C, d) -homogeneous on B for some C. 

We want rates of convergence which hold uniformly over all regions where p is doubling. We 
therefore also require (Definition [3) that C and d from Definition [2] are uniformly upper bounded. 
This will be the case in many situations including the above examples. 

Definition 3. The marginal p is (Co, do) -maximally-homogeneous for some Co > 1 and do > 1, 

if the following holds for all x £ X and r > 0: suppose there exists C > 1 and d > 1 such that p is 
(C, d)-homogeneous on B{x, r), then p is (Co, do)-homogeneous on B{x, r). 

We note that, rather than assuming as in Definition [3] that all local dimensions are at most do, we 
can express our results in terms of the subset of X where local dimensions are at most do- in this 
case do would be allowed to grow with n. The less general assumption of Definition [3] allows for a 
clearer presentation which still captures the local behavior of fc-NN regression. 
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3 Overview of results 



3.1 Local rates for fixed k 

The first result below establishes the rates of convergence for any k > In n in terms of the (unknown) 
complexity on B(x, r) where r is any r satisfying p(B(x, r)) > Q(k/n) (we need at least fi(fc) 
samples in the relevant neighborhoods of x). 

Theorem 1. Suppose p, is (Co, d$)-maximally-homogeneous, and B has finite VC dimension Vb- 
Let < 5 < 1. With probability at least 1 — 25 over the choice of (X, Y), the following holds 
simultaneously for all x G X and k satisfying n > k > Vb In 2n + bx{8/8). 

Pick any x € X. Let r > satisfy p(B(x,r)) > SCak/n. Suppose p is (C, d)-homogeneous on 
B(x, r), with 1 < C < Cq and 1 < d < do. We have 



K(l) k \np(B(x,r)) J 

Note that the above rates hold uniformly over x, k > In n, and any r where p(B(x, r)) > ft(k/n). 
The rate also depends on p,(B(x, r)) and suggests that the best scenario is that where x has a small 
neighborhood of large mass and small dimension d. 



3.2 Minimax rates for a doubling measure 

Our next result shows that the hardness of the regression problem is not tied to a particular choice 
of the metric X or the doubling measure /i. The result relies mainly on the fact that p is doubling on 
X. We however assume that p, has the same expansion constant everywhere and that this constant 
is tight. This does not however make the lower-bound less expressive, as it still tells us which rates 
to expect locally. Thus if [i is (C, (i)-homogeneous near x, we cannot expect a better rate than 
O (n~ 2 /( 2+<i )) (assuming a Lipschitz regression function /). 

Theorem 2. Let p be a doubling measure on a metric space (X , p) of diameter 1, and suppose p 
satisfies, for all x € X, for all r > and < e < 1, 

C x e- d p{B{x,er)) < p(B(x,r)) < C 2 e- d p{B(x,er)), 

where C\, C 2 and d are positive constants independent of x, r, and e. Let y be a subset o/R and 
let A > 0. Define D^x as the class of distributions on X x y such that X ~ p and the output 
Y = f(X) +Af(0, 1) where f is any X-Lipschitz function from X to y. Fix a sample size n > and 
let f n denote any regressor on samples (X, Y) of size n, i.e. f n maps any such sample to a function 
/n|(x,Y) (') ■ X ^ y in L 2 (p). There exists a constant C independent of n and A such that 

. f E X,Y, t |/n|(X,Y)pE) - f{x)\ 2 

{MvTx \2d/(2+d) n -2/(2+d) -°- 



3.3 Choosing k for near-optimal rates at x 

Our last result shows a practical and simple way to choose k locally so as to nearly achieve the 
minimax rate at x, i.e. a rate that depends on the unknown local dimension in a neighborhood 
B(x, r) of x, where again, r satisfies p(B(x, r)) > Vl(k/n) for good choices of k. It turns out that 
we just need p(B(x, r)) > il(n -1 / 3 ). 

As we will see, the choice of k simply consists of monitoring the distances from x to its nearest 
neighbors. The intuition, similar to that of a method for tree-pruning in [7j, is to look for a k that 
balances the variance (roughly 1/fc) and the square bias (roughly n (x)) of the estimate. The 
procedure is as follows: 

Choosing k at x: Pick A > max; p [x, Xi), and pick 8 n s > Inn/ 5. 
Let ki be the highest integer in [n] such that A 2 • 9 n j/ki > n {x). 

Define k 2 = ki + 1 and choose k as argminj.. j S [ 2 ] I n ,s/ki + rf . n (x) ) . 
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The parameter 6 n g guesses how the noise in Y affects the risk. This will soon be clearer. Perfor- 
mance guarantees for the above procedure are given in the following theorem. 
Theorem 3. Suppose /i is (Co, do) -maximally -homogeneous, and B has finite VC dimension Vg. 
Assume k is chosen for each x £ X using die above procedure, and let f n ,k{x) be the corresponding 
estimate. Let < 8 < 1 and suppose n 4/ ( 6+3d o) > (y B In 2n + ln(8/<5)) /9 n ,s- With probability at 
least 1 — 2(5 over the choice o/(X, Y), the following holds simultaneously for all x G X. 

Pick any x € X. Let < r < A satisfy (j,(B(x, r)) > QConT 1 ^. Suppose p, is (C, d)-homogeneous 
on B{x, r), with 1 < C < Co and 1 < d < da. We have 

l /MW . /W |.<(^ + «.) (1+ ^(_^d_) VM , 

where C n>5 = (V B • t^(S/2n) ■ hx(2n/S) + a\) K(0)/K(1). 

Suppose we set 9 ni s = In 2 n/S. Then, as per the discussion in Section |2~T] if the noise in Y is 
Gaussian, we have ^(5/2n) = 0(\nn/S), and therefore the factor C n> s/O n ,& — 0(1). Thus 
ideally we want to set d n ,s to the order of (ty(5/2n) • lnn/S). 

Just as in Theorem [1] the rates of Theorem [3] hold uniformly for all x E X, and all < r < A 
where fi(B(x, r)) > i7(n _1 / 3 ). For any such r, let us call -B(a;, r) an admissible neighborhood. It is 
clear that, as n grows to infinity, w.h.p. any neighborhood B(x, r) of x, < r < sup x , eX p (x, x'), 
becomes admissible. Once a neighborhood B(x, r) is admissible for some n, our procedure nearly 
attains the minimax rates in terms of the local dimension on B(x,r), provided /i is doubling on 
B(x, r). Again, the mass of an admissible neighborhood affects the rate, and the bound in Theorem 
|3]is best for an admissible neighborhood with large mass fj,(B(x, r)) and small dimension d. 



4 Analysis 



Define f n ,k(x) = ^y|x fn,k(%) = J2i w i,k{x)f(Xi). We will bound the error of the estimate at a 
point x in a standard way as 



\fn,k( X ) ~ f( X )\ ^ 2 fn,ki x ) ~ fn,ki x ) 



+ 2 



fn,k(x) - f(x) 



(1) 



Theorem[T]is therefore obtained by combining bounds on the above two r.h.s terms (variance and 
bias). These terms are bounded separately in Lemma|2]and Lemma[3]below. 



4.1 Local rates for fixed k: bias and variance at x 

In this section we bound the bias and variance terms of equation (fl~|i with high probability, uniformly 
over x E X. We will need the following lemma which follows easily from standard VC theory lfl6l 
results. The proof is given as supplement in the appendix. 

Lemma 1. Let B denote the class of balls on X, with VC-dimension Ve- Let < 5 < 1, and define 
Oi n = (Vg In 2n + ln(8/£)) jn. The following holds with probability at least 1 — 8 for all balls in 
B. Pick any a > a n . Then fi(B) > 3a =>• fJ, n (B) > a and ii n (B) > 3a =>■ p(B) > a. 

We start with the bias which is simpler to handle: it is easy to show that the bias of the estimate 
at x depends on the radius rk, n ( x )- This radius can then be bounded, first in expectation using the 
doubling assumption on fi, then by calling on the above lemma to relate this expected bound to 
fk,n( x ) w i m high probability. 

Lemma 2 (Bias). Suppose (i is (Co, dif)-maximally-homogeneous. Let < 6 < 1. With probability 
at least 1 — 5 over the choice o/X, the following holds simultaneously for all x G X and k satisfying 
n > k > V B ln2n + ln(8/<5). 

Pick any x € X. Let r > satisfy /j,(B(x, r)) > 3Cok/n. Suppose /j is (C, d)-homogeneous on 
B(x, r), with 1 < C < Co and 1 < d < do. We have: 
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Proof. First fix X, x £ X and k € [n] . We have: 



/n,fc(») - /(z) 



< 



< }^W it k(x)Xp (Xj,x) < \r k , n (x). 



(2) 



We therefore just need to bound 7** n (a;). We proceed as follows. 

Fix a; € A" and k and pick any r > such that fi(B(x,r)) > 3Cok/n. Suppose /i is (C,d)- 
homogeneous on B(x, r), with 1 < C < Co and 1 < d < do- Define 



3Ck 



1/d 



so that e < 1 by the bound on fi(B(x,r)); then by the local doubling assumption on B(x,r), 
we have fi(B(x,er)) > C~ 1 e d [j,(B(x,r)) > 3k/n. Let a n as defined in LemmaQ] and assume 
k/n> a n (this is exactly the assumption on k in the lemma statement). By LemmaQ] it follows that 
with probability at least 1 — 5 uniform over x, r and k thus chosen, we have /i n ((B(x, er)) > k/n 
implying that r kin (x) < er. We then conclude with the lemma statement by using equation (g). □ 

Lemma 3 (Variance). Let < 5 < 1. With probability at least 1 — 25 over the choice of (X, Y), 
the following then holds simultaneously for all x £ X and k £ [n]: 



fn,k(x) - fn,k(x) 



< 



K(Q) V B -t 2 Y (5/2n)-\n(2n/5) + <T\ 
K{1) k 



Proof. First, condition on X fixed. For any x £ X, k £ [k], let Yj.^ denote the subset of Y 
corresponding to points from X falling in B(x, rk, n {x)). For X fixed, the number of such subsets 
Y x , k is therefore at most the number of ways we can intersect balls in B with the sample X; this is 
in turn upper-bounded by ti Vb as is well-known in VC theory. 

Let ip(Y Xyk ) = fn,k( x ) ~ fn.k(x) ■ We'll proceed by showing that with high probability, for all 
x £ X, ip(Y X! k) is close to its expectation, then we bound this expectation. 

Let 5q < l/2n. We further condition on the event 34 that for all n samples Yi, \Yi — f(Xi)\ < 
ty(5o)- By definition of £y(#o), tne event yg happens with probability at least 1 — n5o > 1/2 . It 
follows that for any x £ X 

E^(Y X)fc ) > ¥(y So ) ■ E ^(Y x>fc ) > I E i>(Y x , k ), 



2y Sn 



where Ey 5 [•] denote conditional expectation under the event. Let e > 0, we in turn have 



1 (3s, k, *l){Y x , k ) > 2E^{Y x , k ) + e) < P ( 3x, k, ^{Y x , k ) > E ^{Y x , k ) + e 



< 



»y So hx,k, ip{Y Xtk ) > E iP(Y x , k ) + e^j +n5 . 

This last probability can be bounded by applying McDiarmid's inequality: changing any Y{ value 
changes tp(Y x . k ) by at most Wi yk ■ ty(<5o) when we condition on the event y$ . This, followed by a 
union-bound yields 

P^ o (3x, k, ij(Y x , k ) > E tp(Y X} k) + < n Vfi cxp j-2e 2 /4(<5 ) £ u£ fc J . 
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Combining with the above we get 



¥(3x € X, ip(Y Xt k) > 2EV(Y x>fe ) + e) < n Vfi exp S^-2e 2 /t Y (6 ) £ w 2 fc | + 
In other words, let <5q = #/2n, with probability at least 1 — 5, for all a; € A" and fc € [n] 



n(5 . 



/n,fc(a:) - /n,fc(a;) 



< 



Y|X 



fn,k( x ) ~ fn,k( x ) 



t 2 Y (S/2n)[V B ln(2n/5)J2^lk 



< 8 E 

Y|X 



fn,k{x)-f n ,k[x) +t 2 Y (S/2n) (v B ln(2n/5)J2wl 



k ! 



where the second inequality is an application of Jensen's. 

We bound the above expectation on the r.h.s. next. In what follows (second equality below) we use 
the fact that for i.i.d random variables Zi with zero mean, E |^ Zi\ = J^i ^ \zi\ ■ We have 



E 

YIX 



fn,k( X ) - fn,k( X ) 



= E 
YIX 



£><,fe(aO (Yi-f(Xi)) 

i 

J2wl k (x)E \Yi - fiX^ 2 < 



Combining with the previous bound we get that, with probability at least 1 — 8, for all x and k, 

U,k(x) - fn, k (x)\ 2 < (V B -%(5/2n)-ln(2n/6)+a Y )-J2 w tk(x)- (3) 

i 

We can now bound J^. w 2 k {x) as follows: 



> w 2 k (x) < m&xwi k(x) = max 



K(p(x, Xi)/r k , n {x)) 



< 



K(0) 



< 



»eW J2jK (p(x,Xj)/r ki n(x)) J2j K (p(x,Xj)/r k , n (x)) 
K(0) 



< 



K(0) 

J2 Xj eB(x,r k:n (x)) K (p( X ^ X j)/ r Kn( X )) ~ K(l)k' 

Plug this back into equation [3] and conclude. 
4.2 Minimax rates for a doubling measure 



□ 



The minimax rates of theorem [2] (proved in the appendix) are obtained as is commonly done by 
constructing a regression problem that reduces to the problem of binary classification (see e.g. [1 
|2][T0)). Intuitively the problem of classification is hard in those instances where labels (say — 1, +1) 
vary wildly over the space X, i.e. close points can have different labels. We make the regression 
problem similarly hard. We will consider a class of candidate regression functions such that each 
function / alternates between positive and negative in neighboring regions (/ is depicted as the 
dashed line below). 



+ 



The reduction relies on the simple observation that for a regressor /„ to approximate the right / 
from data it needs to at least identify the sign of / in the various regions of space. The more we can 
make each such / change between positive and negative, the harder the problem. We are however 
constrained in how much / changes since we also have to ensure that each / is Lipchitz continuous. 
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4.3 Choosing fc for near-optimal rates at x 

Proof of Theorem\3\ Fix x and let r, d, C as defined in the theorem statement. Define 

n ' 5 V 3C J \nfi(B{x,r 
Note that, by our assumptions, 

2/(2+d) 

/i(5(a:, r)) > 6C(9„ ^n" 1 / 3 > 6C(9„ 5 n- d /( 2+d ) = 6C6» n 4 > 6C-. (4) 

71 71 

The above equation implies e < 1. Thus, by the homogeneity assumption on B(x,r), 
fi(B(x, er)) > C~ 1 e d fi(B(x, r)) > 3n/n. Now by the first inequality of (0]i we also have 

1 > ^£ n 4/(6+3d) > ^4/(6+3*) > 
71 71 71 

where a„ = (Va ln2?i + ln(8/<5)) /n is as defined in LemmaQ] We can thus apply Lemma[T]to 
have that, with probability at least 1 — 6, [i n (B(x, er)) > n/n. In other words, for any fc < k, 
fk,n(x) < er. It follows that if k < k, 

Remember that the above inequality is exactly the condition on the choice of fci in the theorem 
statement. Therefore, suppose k% < k, it must be that fc 2 > k otherwise fc 2 is the highest integer 
satisfying the condition, contradicting our choice of k\. Thus we have (i) 6 rh s/k2 < 8 n ,s/ 'k = e 2 . 
We also have (ii) r^ 2n {x) < 2 x l d er. To see this, notice that since k\ < k < fc 2 = fci + 1 we have 
&2 < 2k; now by repeating the sort of argument above, we have fi(B(x, 2 1 / d er)) > 6/t/n which by 
Lemma Q] implies that fi n (B(x, 2 x / d er)) > 2n/n> k^jn. 

Now suppose instead that k\ > k, then by definition of k\, we have (iii) 

/ >9 A 2 ' @n S A 2 ■ 6 n 5 , . , 9 

r kl ,n(x) 2 < < ^ - Ae 2 . 

fci K 

The following holds by (i), (ii), and (iii). Let k be chosen as in the theorem statement. Then, whether 
ki > k or not, it is true that 



(jf <(1 + 4AV = (1 + 4A 2 ) (- 



3C6 nS ^ 2/(2+^^) 



n^x{B(x,r)) 

Now combine Lemma[3]with equation (O and we have that with probability at least 1 — 28 (account- 
ing for all events discussed) 

\U,k{x) - f(x)\ 2 < 2 ~^f d -f + 2A 2 7i„(x) < {^f + 2A 2 ) (°f + rl n (x) 



V J \n^(B(x,r)) 



2/(2+d) 



□ 



5 Final remark 

The problem of choosing fc = k(x) optimally at x is similar to the problem of local bandwidth 
selection for kernel-based methods (see e.g. [17 18 1), and our method for choosing fc might yield 
insights into bandwidth selection, since fc-NN and kernel regression methods only differ in their 
notion of neighborhood of a query x. 
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Appendix 



A Proof of Lemma Q] 

Lemma 4 (Relative VC bounds |[T6l ). Let B be a class of subsets of X. Let < S < 1. Suppose 
a sample X of size n is drawn independently at random from a distribution \x over X with resulting 
empirical distribution /!„. Define a n = (Vg In 2n + ln(8/<5)) /n. 

Then with probability at least 1 — 5 over the choice o/X, all B £ B satisfy 

KB) < fi n (B) + \J )i n (B)a n + a n , and 
Mn(S) < p,(B) + ^n(B)a n + a n . 



Lemma[T]is then obtained as the following corollary to Lemma|4]above. 

Corollary 1. Let B denote the class of balls on X. Let < 6 < 1, and as in Lemma^\above, define 
a„ = (V B ln2n + ln(8/5))/n. 

The following holds with probability at least 1 — 8 for all balls in B. Pick any a > a n . 

• mC- 8 ) > 3a /"n(-B) + y/^n{B)a n > fx(B) - a n > a + y/aa n /i„(B) > a. 

• (J-n(B) > 3a /i(-B) + y/n(B)a n > fx n {B) - a n > a + ^aa n > a. 



B Proof of Theorem |2] 



The minimax rates shown here are obtained as is commonly done by constructing a regression 
problem that reduces to the problem of binary classification (see e.g. 0] |2] [TO]). Intuitively the 
problem of classification is hard in those instances where labels (say — 1, +1) vary wildly over the 
space X, i.e. close points can have different labels. We make the regression problem similarly 
hard. We will consider a class of candidate regression functions such that each function / alternates 
between positive and negative in neighboring regions (/ is depicted as the dashed line below). 



+ 



The reduction relies on the simple observation that for a regressor /„ to approximate the right / 
from data it needs to at least identify the sign of / in the various regions of space. The more we can 
make each such / change between positive and negative, the harder the problem. We are however 
constrained in how much / changes since we also have to ensure that each / is Lipchitz continuous. 
Thus if / is to be positive in some region and negative in another, these regions cannot be too close. 
We therefore have to break up the space X into a set of regions whose centers are far enough to 
ensure smoothness of the candidates /, yet in a way that the set is large so that each / could be 
made to alternate a lot. We will therefore pick the centers of these regions as an r-net of X for some 
appropriate choice of r; by definition the centers would be r far apart and together they would form 
an r-cover over space. As it turns out, any r-net is large under some tightness conditions on the 
expansion constants of \x. 

We start with the following lemma which upper and lower bounds the size of an r-net under some 
tightness conditions on the doubling behavior of /i. Results of this type appear in different forms in 
the literature. In particular similar upper-bounds on r-net size are discussed in |[T3l[T4l . Here, we 
are mainly interested in the lower-bound on r-net size but show both upper and lower bounds for 
completion. Both upper and lower bounds re lie on union-bounds over sets of balls centered on a net. 
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Lemma 5. Let p be a measure on X such that for all x G X, for all r > and < e < 1, 

Cx€- d n{B{x,er)) < p(B(x,r)) < C 2 e- d p{B(x,er)), 

where Ci, C2 and d are positive constants independent of x, r, and e. Then, there exist C[ , C 2 such 
that, for all x G X, for all r > and < e < 1, an (er)-net of B(x, r) has size at least C' 1 e~ d and 
at most C^e . 

Proof. Fix B(x, r), and consider an (er)-net Z of B(x, r). First we handle the upper-bound on \Z\. 

Since, by a triangle inequality, we have for any z G Z, B(x,r) C B(z,2r), it follows by the 
assumption on p that, 



B (z, |r) J > C-H- d ( d ^i (B(z, 2r)) > C^ 1 'i' d e d fi (B(x, r)) 



Now, since for any such z, B (z, |r) C B(x, 2r), and the balls B (z, |r) for z e Z are disjoint 
(centers are er apart) we have 

\Z\ C^4- d e d v (B(x, r)) < /* ( |J B (z, |r) J < /i(B(s, 2r)) < C 2 2 d p(B(x, r)), 

\zez J 

implying that |Z| < C%8 d e- d . 

The lower-bound on \ Z\ is handled similarly as follows. We have for all z G Z 

^(S(z,er)) < C^eVflfor)) < Cf 1 e <, /i(B(a: 1 2r)). 
Thus, applying a union-bound we obtain 

|Z| Cf^^^Ca:, 2r)) > /i [ |J B (z, er) j > p(B{x, r)) > C^^Bix, 2r)), 

implying that |Z| > C^"^. □ 
The proof of the minimax theorem follows. 

Proof of Theorem^ In what follows, for any function g : A" M> M, define ||g|| 2 = K x \g(x)\ 2 , i.e. 
the squared norm of g in L 2 (p). 

Let r„ = (A 2 n) _1 ^ 2+d ^, and let Z be an r„-net of X. Define r = min ||Cy d , ||. For every 
z E Z we define the following function 

A 
5 

It is easy to verify that for all x, x' G X 



f z (x) = - (rr n - p (x, z)) + i.e. f z (x) = whenever p (x, z) > rr n 



\fz{x) - f z {x')\ < ^ \p(x,z) -p(x',z)\ < ^p(x,x'), 
5 5 

that is, f z is ^-Lipschitz. 

Now consider the random vector <r = {s z } zeZ wnere the coordinates <r 2 G { — 1,1} are independent 
Bernoulli r.v.s that are 1 with probability 1 /2. For every possible value of <r define 



M x ) = ^ z g z (x). 



z£Z 



Next we verify that / s is A-Lipschitz. First pick x and x' from the same ball B(z, rr n ). It is clear 
that 
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Now suppose x or x' or both are outside all balls {B(z,rr n )} zeZ , then again it is easy to see 
that \fq(x) — AC^OI < Xp(x,x') . We now check the final case where x € B(z,rr n ) and x' € 
B(z' , rr n ), 2 7^ z'. To this end, first notice that the ring B(z, r„/2) \ B(z, rr n ) is non-empty since 
p{B(z, r n /2)) > Ci(2r)^ d p,(B(z, rr„)) > p,(B(z, rr n )). Pick x" in this ring, and notice that x" 
is outside both balls B(z, rr n ) and B(z' , rr n ). Thus, g z (x") = g' z ( x ") = 0, and we can write 

lAC*) - A(V)I = - <^(V)I < \gM - 9z(x")\ + W) - g'M 

<±{p (x, x") + p (i", a/)) < ^ (2p (z, z") + p (x, x')) 
5 5 

< j (2r„ + p (a, a/)) < ^ (4p (x, x') + p (a, x')) = Xp {x, x') . 
5 5 

At this point we can define T> as the class of distributions on X x y, where X ^ p and Y = 
APO + jV(0, 1), for some A as constructed above. Clearly 2? C 2? m .a and we just have to show 
that 

/n f SUP V2Z2 ^ °(!)- 

Un} v A r„ 

Fix a regressor /„, that is /„ maps any sample (X, Y) to a function in L 2 (p), which, for simplicity 
of notation, we also denote by f n . For (X, Y) fixed, we denote by f n z the projection (in the 
Hilbert space L 2 (p)) of /„ onto the orthonormal system {A/ \\f z \\} z£Z - ^ n °t ner words, f n< z = 
Y. z ez < \\}lf > fz = J2 ze z w n,zf z - Thus, for A fixed, we have 

x E Y || A - All 2 > X E Y \\fn,z - All 2 = x E Y £ - • 

Z<EZ 

To bound ||/ 2 ||, notice that f z is at least Arr„/10 on the ball rr„/2). This ball in turn has mass 
at least C 2 ~ 1 (Tr n /2) d , so ||/ z || > C3/ '-y/n, for C3 appropriately chosen. We therefore have 

C 2 C 2 
sup E ||/„ - / s || 2 > — sup E S~] (w nz — <r z ) 2 > — sup E V 1 {w n , z ■ & < 0} 
X,Y n x, X,Y ^ ri x> X > Y ^rl 

z£Z z£Z 

c 2 c 2 

> — E E Vl {w n . z • & < 0} > -3-E V E 1 {w„, 2 • ^ < 0} . 

For 2 S Z fixed, Ey^x 1 \w n ,z ■ < 0} is the probability of error of a classifier (which outputs 
sign(w„ z )) for the following prediction task. Let xm , X(v\ , • • • %{m) denote the values of X falling 
in B(z, rr n ) where f z is non zero. Then 

(X(l)>- ■■ Y {m)) = S z (fz{X(i)), . . ■ , A(^(m))) + -A/ r (0,i m ) 

is a random vector sampled from the equal-weight mixture of two spherical Gaussians in M. m cen- 
tered at u = (f z (x[i)) 1 . . . , fz(x( m ))) and —u. The prediction task is that of identifying the right 
mixture component from the single sample (Yn\ , . . . Yr m ) ) ■ The smallest possible error for this task 

is that of the Bayes classifier and is well known to be $(— ||m||) > $ (—y/J2i=i / 2 (^«)) ■ Since 

is convex, we can apply Jensen's inequality as follows. 



i=l 



^y $ 

z£Z 



■\lHfz\\ 



The norm ||/ z || is at most C' 3 /y/n since A Arr„ everywhere, and non zero only on the ball 
B(z, rr„) which has mass at most C 1 " 1 (rr„) d . Finally, remember that by Lemma|5] \Z\ > C[r~ d . 
Using these two facts we have 

T A^x*V ll/n - AH 2 > ^ Cf^*(-CS). 

which concludes the proof since /„ is arbitrarily chosen. □ 



12 



